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Abstract. We present the detection of the ammonia (NH3) (J,K) = (1,1) to (4,4) and (6,6) inversion lines toward the prototypi- 
cal Seyfert 2 galaxy NGC 1068, made with the Green Bank Telescope (GBT). This is the first detection of ammonia in a Seyfert 
galaxy. The ortho-to-para-NH^ abundance ratio suggests that the molecule was formed in a warm medium of at least 20 K. For 
the NH3 column density and fractional abundance, we find (1.09±0.23)xl0 14 cnT 2 and (2.9±0.6)xl0~ 8 , respectively, from the 
inner ~1.2 kpc of NGC 1068. The kinetic temperature can be constrained to 80+20 K for the bulk of the molecular gas, while 
some fraction has an even higher temperature of 140±30 K. 
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1. Introduction 

Studies of molecular gas provide information about the gas 
density, the temperature, and kinematics within galaxies, help 
us to understand chemical evolution over cosmic time and al- 
low us to study the triggering and fueling mechanisms of star 
formation and active galactic nuclei (AGN). Most stars are 
formed in dense gas cores, which are embedded in giant molec- 
ular clouds (GMCs). The star-forming activity is related to the 
dense gas and not to the bulk of the GMCs material (Gao 
& Solomon 2004). Determining the physical properties of the 
dense gas in galaxies is therefore of fundamental importance 
for our understanding of star formation and the evolution of 
galaxies. Among the most commonly observed species are CO, 
CS, HCN, and HCO + . In local dark clouds, the temperature can 
be constrained by observations of the 7=1-0 transition of CO, 
both because this transition is opaque and thermalized and be- 
cause the emission is often extended enough to fill the beam of 
a single-dish telescope. However, in external galaxies, filling 
factors are much less than unity. Furthermore, CO, CS, HCN, 
and HCO + suffer from a coupled sensitivity to the kinetic tem- 
perature and spatial density, making an observed line ratio con- 
sistent with both a high density at a low temperature and a low 
density at a high temperature. Specific information about the 
individual physical parameters therefore requires a molecular 
tracer that possesses a singular sensitivity to either density or 
temperature. Ammonia (NH3) is such a molecule. 

Ammonia is widespread and its level scheme contains in- 
version doublets owing to the tunneling of the nitrogen atom 
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through the plane defined by the three hydrogen atoms. The 
metastable ( J, K — J) rotational levels, for which the total an- 
gular momentum quantum number J is equal to its projection 
on the molecule's symmetry axis, are not radiatively but col- 
lisionally coupled. NH3 is therefore a good tracer of the gas 
kinetic temperature. Another advantage is that the inversion 
lines are quite close in frequency, thus allowing us to measure 
sources with similar beam sizes and to use the same telescope- 
receiver combination, which minimizes calibration uncertain- 
ties of the measured line ratios. As a consequence, NH3 is 
widely studied to investigate the physical properties of dark 
clouds and massive star-forming regions in our Galaxy (e.g., 
Ho & Townes 1983; Walmsley & Ungerechts 1983; Bourke et 
al. 1995; Ceccarelli et al. 2002; Pillai et al. 2006; Wilson et al. 
2006). However, in spite of the high sensitivity and stable base- 
lines of present state-of-the-art facilities, ammonia multilevel 
studies in extragalactic sources are still rare and are limited to 
the Large Magellanic Cloud, IC 342, NGC 253, M 51, M 82, 
Maffei 2, and Arp 220 (Martin & Ho 1986; Henkel et al. 2000; 
WeiB et al. 2001; Takano et al. 2002; Mauersberger et al. 2003; 
Ott et al. 2005, 2010; Takano et al. 2005 ) in the local Universe, 
and in absorption to the gravitationally lensed B0218+357 and 
PKS 1830-211 at redshifts of z ~ 0.7 and 0.9 (Henkel et al. 
2005, 2008). 

NGC 1068, a prototypical Seyfert 2 galaxy (Antonucci & 
Miller 1985), is located at a distance of 15.5 Mpc (a heliocen- 
tric systemic velocity cz — 1 137 km s is used throughout this 
paper, adopted from the NASA/IPA Extragalactic Database), 
making it one of the nearest Seyfert galaxies and thus an ideal 
target to investigate the physical properties of the molecular 
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gas in the vicinity of an AGN. Here we report the first detec- 
tion of ammonia in this prominent Seyfert 2 galaxy to evaluate 
the kinetic temperature of its dense molecular gas. 

2. Observations 

We observed ammonia inversion lines toward the nucleus of 
NGC 1068 with the Green Bank Telescope (GBT) of the 
National Radio Astronomy ObservatorjQ on 2005 October 14 
and October 19. We configured the telescope to observe both 
circular polarizations simultaneously using a total-power nod- 
ding mode in which the target position was placed alternately 
in one of the two beams of the 22-26 GHz K-band receiver. 
The full-width-at-half-maximum (FWHM) beam size was ap- 
proximately 31" in the observed frequency range, between 
23.6 and 25.0 GHz. The pointing accuracy was ~ 5". To ob- 
serve the NH3(1,1) through (4,4) lines, we configured the spec- 
trometer with an 800 MHz spectral window centered on the 
NH3(3,3) line. The channel spacing was 390 kHz, correspond- 
ing to ~ 5 kms~'. We averaged the data from the two dates to 
produce the final spectrum, which has an on-source integration 
time of 4 hours. We observed the NH3(6,6) line only on the 
second date, 2005 October 19. Here we used a 200 MHz spec- 
tral window observed with a channel spacing of 24 kHz, and 
then averaged channels in post-processing to achieve 390 kHz 
channel spacing. The total integration time for the NH3(6,6) 
line was 1 hour. 

The data were calibrated and averaged in GBTIDL. For 
calibration, we used an estimate of the atmospheric opacity 
obtained from a weather model. We subtracted a polynomial 
baseline fitted to the line-free channels adjacent to each ammo- 
nia line. For the (1,1) to (4,4) lines observed simultaneously, 
the relative calibration is very good and is limited primarily by 
the uncertainty in the baseline fits, which is 10% for the (1,1) 
to (3,3) lines, 15% for the (4,4) line, and 27% for the (6,6) line, 
respectively. The (1,1) and (2,2) line profiles are close to each 
other but are marginally separated and the Gaussian fit for each 
line does not suffer from this effect. With the calibration error 
of 15% for the GBT telescope itself, we estimate the absolute 
calibration uncertainties, including the error from the baseline 
fits, to be ~ 18-31% for the ammonia lines. 

3. Results 

3.1. NH 3 lines 

This is the first time that ammonia is detected in NGC 1068, 
making it the first Seyfert galaxy observed in this molecule. 
Also, because Arp 220, B0218+357, and PKS 1830-211 were 
detected in absorption, NGC 1068 is so far the most distant 
galaxy where ammonia is detected in emission. Three para, 
(/, K) = (1,1), (2,2), and (4,4), and two ortho, (J, K) = (3,3) and 
(6,6), transitions are covered by the observations, and all five 
are detected. Figure Q] shows the spectra before and after sub- 
tracting a polynomial baseline to the line free channels. The 
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line parameters are listed in Table [T] where integrated flux, /, 
peak flux density, S , central velocity, V, and line width, AV1/2, 
are obtained from the results of Gaussian fitting. 

The line profiles are similar for the (1,1), (2,2) and 
(3,3) lines, which share a similar central velocity of around 
-35 kms" 1 , relative to a heliocentric systemic velocity of 
cz= 1137 km s _1 , and a FWHM line width of 220-274 kms" 1 . 
The similar central velocities and the comparable widths of the 

(1.1) , (2,2) and (3,3) lines suggest that the line emission comes 
from the same region. The (4,4) and (6,6) lines have different 
central velocities and a narrower line width. The (4,4) line is 
located at a frequency where the baseline is particularly steep. 
The (6,6) line is weak and falls into a window with a slightly 
distorted baseline, making it difficult to determine accurate line 
parameters. The comparatively narrow line widths of the (4,4) 
and (6,6) lines may suggest that their emission originates from 
a less extended region than that of the lower excitation lines. 

Our measured continuum flux density of NGC 1068 ranges 
from 0.33 to 0.37 Jy at frequencies of 23.4 to 25.1 GHz. This 
is consistent with the value of 0.342+0.034 at 22 GHz reported 
by Ricci et al. (2006). 

3.2. W/-/ 3 column density and rotation temperature 

Assuming that the line emission is optically thin and the con- 
tribution from the cosmic background is negligible, the sum 
of the beam-averaged column densities of the two states of an 
inversion doublet can be calculated using 

, s 1.55 x 10 I4 /(/+ 1) r 
N(J,K) = J r mb dv (1) 

(e.g., Mauersberger et al. 2003), where the column density 
N(J,K), the frequency v, and the integrated line intensity 
J r m bdv, based on the main beam brightness temperature, r m b, 
are in units of cm 2 , GHz, and K kms 1 , respectively. The cal- 
culated column densities (1 Jy corresponds to 2. 16 K on a r m b 
scal^l) are given in Table [T] 

Following the analysis described by Henkel et al. (2000), 
the rotational temperature (T ro t) between different energy lev- 
els can be determined from the slope, a, of a linear fit in the ro- 
tation diagram (i.e., Boltzmann plot, which is normalized col- 
umn density against energy above the ground state expressed 
in E/k) by T mt = -loge/a s* -0.434/a. Figure [2] shows the 
rotation diagram including the five measured metastable NH3 
inversion lines. Ignoring differences in line shape, we obtain 
7^ = 92^jg K by fitting the para-NH3 species, i.e. the (1,1), 

(2.2) and (4,4) transitions, and T ml = 125^' K for the (3,3) 
and (6,6) ortho-NH3 transitions. As discussed in § 13.11 the 
(1,1), (2,2) and (3,3) line emission may originate in a region 
that is different from that of the (4,4) and (6,6) lines. The for- 
mer three lines can be fitted by T mt = 61^ 2 K, and the latter 
two by T mt = 1 1 K. The rotation temperature between the 
lowest inversion doublets of para-ammonia, (1,1) and (2,2), is 
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Table 1. NH3 line parameters 



(J,K) 


I" 


5" 


V" 


AVi /2 a 


N(J, K) b 




(Jy km s~') 


(mJy) 


(kms- 1 ) 


(kms" 1 ) 


(10 13 cm" 2 ) 


(1,1) 


0.89( 0.07) 


3.26( 0.13) 


-30(6) 


274(17) 


2.52(0.24) 


(2,2) 


0.79( 0.06) 


3.10( 0.14) 


-37(6) 


254(19) 


1.67(0.19) 


(3,3) 


1.19( 0.10) 


5.40( 0.14) 


-37(3) 


220( 8) 


2.22(0.13) 


(4,4) 


0.55( 0.05) 


3.72( 0.14) 


-2(3) 


149( 7) 


0.96(0.07) 


(6,6) 


0.28( 0.04) 


1.88( 0.23) 


-56(9) 


147(27) 


0.43(0.12) 


(5,5) 










0.53 


(0,0) 










2.55 












10.9(2.3) 



The values are derived from Gaussian fits to the spectra. Given velocities are relative to a heliocentric velocity of cz = 1 137 km s" 1 . 
b The column density for the (5,5) state is extrapolated from the (2,2) and (4,4) states using r rot = 1 19 K, while the value for the (0,0) state 
is extrapolated from the (3,3) state by adopting T m[ = 44 K. The total column density, iVtotai, includes the populations of the levels from 
(0,0) to (6,6). Because relative calibration is excellent up to at least the (4,4) transition, the errors given in parenthesis refer exclusively to 
the Gaussian fits. For the cumulative column density, iVtotab however, the given error accounts for the absolute calibration uncertainty (§|2j. 



T mt =44™ K. Without the (1,1) state, the other two para-NH 3 
transitions, (2,2) and (4,4), yield T mt = 119^ K. 

3.3. /VH 3 abundance 

To estimate the populations of the (J, K) = (0,0) and (5,5) states, 
we assume the rotation temperature T rol = 44 K, derived from 
the two lowest inversion doublets, for the (0,0) ground state, 
and T m = 1 19 K from the fit to the (2,2) and (4,4) lines, for the 
(5,5) state. The derived column densities are 2.55xl0 13 cm -2 
for the (0,0) state, which is not a doublet, and 0.53xl0 13 cm -2 
for the (5,5) state, respectively. The former is extrapolated from 
the (3,3) state and the latter from the (2,2) and (4,4) states. The 
resulting total column density of ammonia is ( 1 .09+0.23)x 10 14 
cm" 2 if only considering the populations in the metastable 
states up to (J, K) = (6,6). This yields an ammonia gas mass 
of 68+14 M within a beam size of 31", which corresponds to 
a physical size of 2.3 kpc. To estimate the mass of molecular 
gas in NGC 1068, we convolved the 17" resolution map in CO 
J = 1-0 by Kaneko et al. (1989) to the GBT's resolution of 31" 
and derived an integrated intensity of 103 K kms 1 . Adopting 
a conversion factor of 0.8 M Q (Kkms -1 pc 2 ) -1 , which de- 
scribes ultraluminous galaxies (Downes & Solomon 1998) and 
less conspicuous nuclear starbursts (Mauersberger et al. 2003), 
we obtain a molecular gas mass of 3.5xl0 8 M . Therefore, 
the ammonia abundance relative to molecular hydrogen is esti- 
mated to be (2.9+0.6)xl0~ 8 within a radius of ~1.2 kpc around 
the center of NGC 1068. 

4. Discussion 

The ammonia abundance in NGC 1068 is consistent with 
the values of (1.3-2.0)xl0~ 8 reported by Mauersberger et al. 
(2003) for the central regions of several nearby galaxies such 
as NGC 253, Maffei 2 and IC 342. It is lower than that of a few 
xlO -7 for the central ~ 500 pc of the Milky Way (Rodriguez- 
Fernandez et al. 2001), but higher than the extremely low value 
of 5xlO- 10 determined for M 82 by WeiG et al. (2001). M 82 
is a special case, because it contains a starburst in a late stage 



of its evolution leaving NH3, which is a molecule particularly 
sensitive to UV radiation (e.g., Suto & Lee 1983), in only a 
few molecular cores that are still sufficiently shielded against 
the intense radiation field. 

With the three para- and two ortho-NH3 lines observed, we 
can estimate the ortho-to-para abundance ratio to be R ~ 0.9 
(we do not give an error here because of the extrapolated (0,0) 
column density) as an upper limit because the (0,0) column 
density might be overestimated by extrapolating from the (3,3) 
column density, which could be affected by maser activity (e.g., 
Walmsley & Ungerechts 1983; Ott et al. 2005). A value near or 
below unity is expected, if ammonia was formed in a warm 
medium at a gas kinetic temperature of at least 20 K and only R 
values above ~1.5 would hint at a cool formation temperature 
of ~ 20 K (e.g., Takano et al. 2002). 

The rotation temperatures from the multilevel study of the 
ammonia inversion lines are a good approximation to the ki- 
netic gas temperature only for low (T^ n < 20 K) temperatures 
(Walmsley & Ungerechts 1983; Danby et al. 1988). At higher 
temperatures, the rotation temperature provides a robust lower 
limit to the kinetic gas temperature as long as saturation effects 
do not play a role. 

To determine the kinetic temperature itself, we here use a 
one-component large velocity gradient (LVG) analysis adopt- 
ing a spherical cloud geometry as described by Ott et al. (2005). 
Dahmen et al. (1998) estimated that the velocity gradient 
ranges from 3 to 6 kms 1 pc -1 for Galactic center clouds and 
Meier et al. (2008) found a typical value of 1 to 2 kms -1 pc -1 
for GMCs. Here we adopt a median value of 3 km s pc -1 for 
the LVG models presented in this paper. The ammonia abun- 
dance derived in § I3.3l is adopted, which yields an NH3 abun- 
dance per velocity gradient of [NH3]/(dv/dr) ~ 10~ 8 pc (km 
s- 1 )- 1 ([NH 3 ] =N(NH 3 )/N(H 2 ) is the fractional abundance of 
ammonia). Changing the velocity gradient by a factor of 3 
will affect the derived kinetic temperature by less than 10 K. 
Figure [3] shows NH3 inversion line ratios as a function of the 
H2 density («h,) and kinetic temperature (7\i n ). As long as the 
NH3 lines are optically thin, the ratios are almost independent 
of the gas density and are therefore a good indicator for gas 
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Fig. 1. Metastable ammonia inversion lines observed to- 
ward NGC 1068 (a 2 ooo=02 h 42 ra 40.7 s , £ 2 ooo = -00°00'48"). 
Vertical lines mark a heliocentric systemic velocity of 
cz — 1 137 km s _1 . Velocities are displayed relative to this value 
(in the upper left panel, the velocity scale refers to the (1,1) 
line). Left: Spectra without baseline subtraction, overlaid with 
a polynomial baseline fit to the line-free channels. Right: Final 
baseline-subtracted spectra overlaid with Gaussian fits. 



kinetic temperature. At high densities the (1,1) and (2,2) lines 
saturate, which causes the steeply declining T^\ n values with 
rising density for a given (2,2)/(l,l) line ratio. While the op- 
tical depth may be well above unity in molecular cores (see, 
Giisten et al. 1981 for the Galactic center region), the bulk 
of the ammonia emission should be optically thin and arises 
from gas densities «h, < 10 5 cm -3 (e.g., Tieftmnk et al. 1998). 
For a typical gas density, tiu 2 , of 10 30 ~ 4 - 8 crrT 3 , the observed 
(2,2)/(l,l) line ratio constrains the kinetic temperature 7\i n to 
80+20 K, which characterizes the bulk of the molecular gas. 
The (6,6)/(4,4) line ratio yields a gas temperature of 140+30 K 
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Fig. 2. Rotation diagram of metastable ammonia transitions to- 
ward NGC 1068 (see § 13.21 1. The open squares show the nor- 
malized column densities determined from the integrated line 
intensities. The numbers mark the rotational temperatures in 
K. The absolute calibration uncertainties, including the domi- 
nant contribution from the baseline fits as well as uncertainties 
in the Gaussian fits and in the overall calibration uncertainty, 
have been taken as error bars (§|2]i. Note, however, that relative 
calibration is excellent up to at least the (4,4) transition. 



for a gas density within the range of 10 3 - 0-5 - 5 cirT 3 . The lat- 
ter indicates the existence of a hotter component of the gas, as 
is also suggested by the high rotational temperature revealed 
by the rotation diagram (Fig. [2). Here the (6,6)/(3,3) line ratio 
is not used to estimate the temperature because the (3,3) line 
emission is partly affected by the maser activity in some pa- 
rameter ranges (Walmsley & Ungerechts 1983), and because 
the (6,6) line was not observed simultaneously as the (3,3) line 
and the uncertainty of the baseline fit of the (6,6) line is largest 
among all spectra as described in §|2] 

The GBT beam size covers the two inner spiral arms 
and the circumnuclear disk (CND) of NGC 1068 (see, e.g., 
Schinnerer et al. 2000). Our observed ammonia line emission 
likely has contributions from the spiral arms and the CND. The 
CND may be warmer because it is heated not only by young 
stars but also by the AGN. There are observations that directly 
support higher temperatures in the CND. Infrared rotational 
transitions of molecular hydrogen (Lutz et al. 1997) indicate 
a wide range of gas temperatures between 100 and 800 K. The 
dust near the AGN has been particularly thoroughly studied 
(e.g., Tomono et al. 2006; Poncelet et al. 2007; Raban et al. 
2009). Because densities may be high within the central arcsec- 
ond 10 5 cm -3 ), the dust and the gas phase may be coupled 
and thus these results may be relevant for the kinetic temper- 
ature of the gas component. With mid-infrared (MIR) multi- 
filter data, Tomono et al. (2006) obtained a dust temperature of 
-200 K within a 1.0"-sized region (-80 pc). 
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Fig. 3. NH3 inversion line ratios as a function of the H2 den- 
sity («h,) and kinetic temperature (T^ n ). The red and blue lines 
mark the (2,2)/(l,l) and (6,6)/(4,4) line ratios, respectively, ob- 
tained from the large velocity gradient model outlined in § [4] 
Dashed lines show the line ratios measured with the GBT. 
[NH 3 ]/(dv/dr) = 1(T 8 pc (km s' 1 )' 1 ([NH 3 ] = N(NH 3 )/N(H 2 ) is 
the fractional abundance of ammonia). At high densities the 
(1,1) and (2,2) lines saturate, which causes the steeply declin- 
ing Tkin values with rising density for a given (2,2)/(l,l) line 
ratio. 

Our data provide a direct estimate of the gas properties for 
molecular gas in the inner 1.2 kpc of the galaxy. However, 
higher angular resolution data are needed to separate the dif- 
ferent components and to isolate the CND in this galaxy. 

5. Conclusions 

Our main results from the ammonia observations toward the 
prototypical Seyfert 2 galaxy NGC 1068 are: 

(1) The metastable NH 3 (J,K) = (1,1) to (4,4) and (6,6) in- 
version lines are for the first time detected in emission toward 
NGC 1068. This opens up a new avenue to determine kinetic 
temperatures of the dense gas in nearby highly obscured AGN. 

(2) For the NH 3 column density and fractional abundance, 
we find (1.09+0.23)xl0 14 cirT 2 and (2.9+0.6)xl0~ 8 in the in- 
ner -1.2 kpc of NGC 1068. 

(3) With an ortho-to-para-NH 3 abundance ratio of ~0.9, the 
ammonia should have been formed in a warm medium of at 
least 20 K. 

(4) The kinetic temperature can be constrained to 80+20 K 
for the bulk of the molecular gas, while some gas fraction has 
an even higher temperature of 140+30 K. 
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